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Introduction: In many applications of perturbative field theory at finite 
temperature it is necessary to compute higher-order corrections to Green 
functions, i.e., loop diagrams. In one approach JIJ, the real-time formal- 
ism, one computes Green functions in real time, i.e., thermal averages of 
time-ordered products of Heisenberg field operators. The perturbative eval- 
uation of these Green functions entails a duplication of the degrees of free- 
dom (through the introduction of ghost fields), which either increases the 
multiplicity of diagrams or forces the use of matrix-valued propagators. In 
the Euclidean or imaginary-time formalism, on the other hand, one considers 
the so-called temperature Green functions, i.e., thermal averages of r-ordered 
Heisenberg field operators in imaginary time t = —ir, with < r < 1/T. 
The Fourier transforms of these temperature Green functions are defined only 
for discrete values of the (Euclidean) energy, P° = u n (n = 0, ±1, . . .), with 
uj n = (2n)irT for bosons and uo n = (2n + 1)itT for fermions, the so-called 
Matsubara frequencies. In this formalism there is no need for ghost fields 
and the Feynman rules are very similar in form to those of zero-temperature 
Euclidean field theory, except for the fact that now each loop in the diagram 
is associated with an infinite sum over Matsubara frequencies, instead of an 
integration over a continuous energy variable. Although these sums can be 
computed in a number of ways, usually in a systematic fashion, these compu- 
tations can become quite tedious for higher loop diagrams. In one procedure, 
the so-called Saclay method 0, one uses a Fourier integral representation of 
the Euclidean propagator in terms of a variable r conjugate to the discrete 
Matsubara loop energy K°. In this form the energy sums are easily done, 
giving rise to delta functions in the r-variables, which, in a graph with / 
internal lines and L loops, can be used to eliminate L of the / r-integrals. 
What remains are r-integrals of exponential functions, each of them giving 
rise to a denominator linear in energies. Although straightforward, in prac- 
tice the calculations usually become cumbersome, and a simple form for the 
final result can only be obtained after due use of some of the identities satis- 
fied by the Bose-Einstein and Fermi-Dirac thermal occupation numbers 
Another way to perform the Matsubara sums is the standard contour inte- 
gration method, in which the summation over the energy is represented as a 
contour integral with an appropiately chosen integrand. In reference this 
method was used, along with a decomposition of the Euclidean propagator 
into partial fractions, to compute general iV-point graphs. Certain rules are 
given that allow to write down the general answer directly in terms of the 
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external energies. In both methods, the result for a general amputated scalar 
graph with N + 1 vertices (N > 1), I internal lines, and external 4-momenta 
Pi = (Ff, Pi), has the form 

^4— / in T^br n (2-) 3 ^ (3) (k y )] d(p°, (i) 

^ J i=l (.27TJ y =1 

where (7 represents the coupling constant, S is the symmetry factor of the 
graph, the E^ = (k| + m?) 1 / 2 are the energies of the internal lines (carrying 
spatial momenta kj), and ky denotes the total momentum entering vertex 
V. The delta functions ensure conservation of spatial momentum at each 
vertex, so that the integration measure reduces essentially to an integration 
over the 3-momenta of the L = I — N independent loops. The integrand 
in (JJ) is of the form Y^aN a /R a , a sum of rational expressions, where each 
numerator N a is a polynomial of degree L in the occupation numbers rij = 
n(Ei) = (exp(_Ej/T) — and each denominator R a is a product of exactly 
N factors, all linear in the energies iP® and Ei. 

In this letter we present simple diagrammatic rules to write down an 
explicit expression for the integrand D(P^,Ei), appearing in ([[]), for any 
finite-temperature Euclidean Feynman graph, thereby completely avoiding 
the need to perform Matsubara sums. Although similar rules have been 
given in Ref. [|J, we believe that ours are simpler and explicitly incorporate 
the known fact that finite temperature amplitudes can be related to zero- 
temperature forward scattering amplitudes in a thermal bath [||]. The rules 
can be formulated in two alternative ways, in terms of two different types 
of diagrams. In one case, the diagrams are closely related to those encoun- 
tered in old-fashioned zero-temperature time-ordered perturbation theory 0, 
and lead to a result that strongly resembles in form the one obtained by us- 
ing standard Euclidean Feynman rules, once the Matsubara sums have been 
performed. A second, simpler form of the rules can be obtained from the pre- 
ceding one by regrouping similar diagrams into standard Feynman diagrams, 
thus reducing drastically the number of diagrams to be considered. In both 
approaches, the diagrams include cut (or "snipped" ) lines that represent the 
interaction of the external particles with the thermal bath. 

We shall first consider diagrams for scalar particles. In this case the 
integrand D depends on the spatial momenta only through the energies of 
the internal lines E^. We will indicate at the end the necessary modifications 
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for fermions. 

The Rules: In our first approach, which is analogous to the old-fashioned 
time-ordered diagrammatic expansion, the Feynman graph for the case at 
hand is given by expression (|l|), with the integrand D(Pf,Ei) computed 
according to the following rules: 

a. For each external line, characterized by a real Euclidean 4- vector (P z °, p/), 
define the energy of the line as zP°. For each internal line define its 
energy as Ei = (k? + m 2 ) 1 / 2 , where k« is the 3-momentum carried by 
the line and m ; is the mass of the propagating particle. 

b. Define a direction of time or energy flow (which we shall take conven- 
tionally from left to right) and consider all possible orderings of the 
vertices along this direction (see, e.g., Figs, l.a and l.b. For a graph 
with N + 1 vertices there will be (N + 1)! such orderings). 

c. For each time-ordered graph generated in (b) consider, in addition to 
itself, all possible connected graphs that can be obtained by snipping 
any number of internal lines. Each line that is snipped becomes a 
pair of legs we shall call thermal legs. Attach a cross to their ends to 
distinguish them from the original external lines of the graph. Both legs 
of a given pair inherit the energy Ei of the internal line that originated 
them. However, one leg must be oriented as incoming with energy Ei 
and the other as outgoing with energy E^. Both possible orientations 
have to be considered, each one generating a different diagram (see, 
e.g., Figs, l.c and l.d). 

d. For each graph in (c), define its total incoming energy, Ei nc , as the 
sum of all incoming external energies plus the energies of all incom- 
ing thermal legs that join the diagram before their outgoing partner 
(e.g. Figs, l.d and l.f). Thermal leg pairs that satisfy this property 
shall be referred to as external, and those that do not as internal 
(e.g. Figs, l.c and l.e). Then, associate to this graph an expression 
equal to the product of the following factors: 

1. draw a full vertical division (a "cut") between each pair of con- 
secutive time-ordered vertices (there are N such cuts in a graph 
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with N + 1 vertices); for each cut include a factor 



where E cut is the total energy of the intermediate state associated 
with the cut, defined as the sum of the energies of all the lines 
that cross the cut in question (as in zero-temperature time-ordered 
perturbation theory), plus the energies of all internal thermal pairs 
whose originating internal line would have crossed the cut. 

2. include a thermal occupation factor ni = n(Ei) for each thermal 
pair (of energy Ei) in the diagram (if any). 

3. include an overall factor of (—1)^, where iV + 1 is the number of 
vertices. 



The integrand D(P®,Ei) is the sum of the expressions computed according 
to rule (d), over all the graphs in (c). 

Since the temperature T enters only through the thermal occupation fac- 
tors Hi described in rule d.2, it is evident that the full result separates explic- 
itly into a zero-temperature and a temperature-dependent part, the former 
being obtained solely from the "un-snipped" diagrams, since the thermal 
occupation factors vanish for T = 0. This is a known property of the Eu- 
clidean finite-temperature result, once the Matsubara sums have been dealt 
with. Rule (c), in conjunction with rule (d.2), embodies the known fact that 
any term contributing to the temperature-dependent part of a L-loop graph 
contains at most L factors of the thermal occupation numbers. In fact, in 
a L-loop graph one can snip at most L internal propagators if the resulting 
graph is to remain connected. 

As it stands, the diagrammatic method gets out of hand for diagrams with 
more than a few vertices, due to the factorial growth in the number of time- 
ordered diagrams (see rule (b)). We show below that the time-ordered graphs 
can be regrouped into standard "covariant" Feynman diagrams, thereby sim- 
plifying the computation considerably. These diagrams are similar to stan- 
dard T = Feynman diagrams, except for the inclusion of thermal leg pairs. 

A general proof of the above diagrammatic rules will be presented else- 
where JFJ. We have verified their correctness in a few simple cases by hand 
and in various nontrivial cases with Mathematica. In the remainder of this 
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letter we would like just to illustrate how the rules work in a few examples. 
We shall leave the analysis of the general case for Ref. [0]. 

Example 1: We consider first the lowest order correction to the inverse 
propagator in the scalar g<f) 3 theory (see Fig. 1), which we will call L 2 (P), 
with P = (P°, p). A straightforward application of the Saclay method gives 



D 2 (P°,E 1 ,E 2 



1 + rii + n 2 



+ 



iP°-E 1 -E 2 iP G -E 1 + E 2 



rii - n 2 

~ iPO + E l -E 2 
1 +n 1 + n 2 

iP° + E 1 + E 2 ' 



with 



where E\ = (k 2 + m 2 ) 1//2 and E 2 = [(p — k) 2 + m 2 ] 1//2 are the energies of 
the internal lines. 

In our approach, each one of the ten terms in equation (|]) corresponds 
to one particular snipped, time-ordered diagram^. The zero-temperature 
contribution is obtained from diagrams l.a and l.b. These give 

£> 2 (P ,-Ei,-E 2 ) = n + ^ (5) 

2V ' u 21 t=o iP°-E 1 -E 2 iP° + E 1 + E 2 KJ 



One can readily verify that the usual T = Feynman integral of this 1-loop 
correction can be put in the general form ([]]), with D 2 given in the equation 
above, by performing explicitly the fco-integral by contour integration. In the 
time-ordered formalism, the relative minus sign between both terms in (^) 
has its origin in the energy of the intermediate state in diagram (2b), namely 
E cut = 2iP° + E x + E 2 . 

The finite temperature piece of D 2 , according to our rules, is obtained 
from the graphs l.c to l.f (and the corresponding ones obtained by inter- 
changing the labels of lines 1 and 2). Their contribution to D 2 is computed 

2 This is true in this particular simple case. In general, there is no such simple corre- 
spondence between the Saclay method and our diagrams. 
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according to rule (d), giving: 



-m 



1 1 

iP° -E l -E 2 + iP° + E X -E 2 



1 

iP° + E 1 + E 2 
1 



iP° ~E 1 + E 2 



(6) 



which fully agrees with the Saclay result (f|). The diagrams l.c to l.f have 
an obvious physical interpretation as extra, thermal contributions to the 
T = diagrams l.a and l.b, where one (or, in general, more than one) of the 
particles participating in the process is taken out from and put back into the 
thermal bath, with probability n(E). 

Example 2: As a less trivial example, consider the 2-loop graph in the A</> 4 
theory, shown in Fig. 2. a. As in the previous case, there are two possible 
time-orderings of the vertices, depicted in Figs. 2. a and 2.b. The T = 
contribution to the integrand D is obtained with our rules directly from 
these graphs, and is equal to: 



1 



+ 



1 



ipo - El -E 2 -E 3 lP o + El + E2 + E . 



(7) 



The finite temperature contribution is obtained by snipping internal lines of 
graphs 2. a and 2.b. As in the previous example, one has terms that are linear 
in the rij's, corresponding to cutting only one internal line. Two particular 
contributions are shown in Figs. 2.c and 2.d. Graph 2.c is calculated with 
E mc = iP° + E l and E cut = E 2 + E 3 to give -rn/(iP° + E 1 - E 2 - E 3 ). 
Similarly, graph 2.d has E inc = iP° and E cut = 2iP° + E\ + E 2 + E 3 , and 
gives +rii/(iP° + E\ + E 2 + E 3 ). A new type of contribution appears in this 
2-loop graph: terms bilinear in the n^'s, obtained by cutting two internal 
lines. Again, two particular contributions are shown in Figs. 2.e and 2.f. For 
graph 2.e, E inc = iP° + E\ + E 3 and E cut = E 2 , which gives —riin 3 / (iP° + 
E 1 -E 2 + E 3 ). For graph 2.f, E mc = iP° + E 3 and E cut = 2iP° + E 1 + E 2 , 
which gives +riin 3 / (iP° + E\ + E 2 — E 3 ). The full finite-temperature result 
is obtained from the sum of all possible snipped diagrams. 

Example 3: As a last example we consider the one-loop contribution to 
the 3-point function in the g<j) 3 theory, shown in Fig. 3. The new ingredient 
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here is the presence of many time-orderings and more than one intermediate 
state. There are 6 possible vertex orderings. For the T = part, two of these 
orderings are shown in Figs. 3. a and 3.b, and their contribution to D is: 

1 



(iP2 -Ek- E 2 )(zP« + iP° -Eh.- E 3 ) 

1 

+ 



zP° + zP 6 ° + E l + E 3 ) (iPO + E l + E 2 ) w 

The finite temperature part comes from snipping the T = graphs. From 
graphs 3. a and 3.b we thus obtain, amongst others, graphs 3.c and 3.d, whose 
contribution to D is: 

n- 2 



(iPS -E 1 + E 2 )(iP* + zP° -E x - E 3 ) 

+ ^ (9) 

^ (zP° + iP b ° + E l + E 3 ) (iP a o + E 1 + E 2 ) { ' 

Again, to obtain the full result one has to consider all possible orderings 
of the vertices and all possible snips of internal lines that keep the graph 
connected. 

In the first example given above, it is easy to realize that the graphs 
in Figs, l.c to l.f can be combined in pairs to form "thermally snipped" 
Feynman diagrams. For instance, the first and fourth terms inside the square 
bracket in @ can be combined into the expression 

2E V°-A) 2 -£f (10) 

which, apart from the factor 2E 2 and the overall minus sign in (H), is precisely 
the Euclidean propagator in the t-channel diagram for the zero-temperature 
forward scattering process 4>{P°, p) + </>(—iEi, k) — > (f>(P°, p) + (p(—iEi, k), 
shown in Fig. 4. a. The factor 2E 2 cancels out a similar factor in the integra- 
tion measure in (P, so that the contributions from graphs l.c and l.f add up 
to 

^TiP^-iE,), (11) 
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where T(P°, —iEi) denotes the t-channel diagram in Fig. 4. a. Denoting like- 
wise the s-channel diagram in Fig. 4.b by S(P°, —iE\), we see that L2(P°, p) 
can be expressed purely in terms of standard Feynman diagrams as: 



L 2 (P°,p) = L 2 (P°,p 
d 3 k 



T=0 ^ 



/ 2E Hl [ T(jP0 ' ~ lEl) + ~ lEl) \ + (1 ~ 2) - (12) 



2#i 

It is not difficult to realize that the ability to regroup a set of our time- 
ordered diagrams into a standard Feynman diagram is quite general. Con- 
sider the sum of all possible time orderings of a given diagram, with a fixed set 
of oriented thermal legs. Alternatively, consider the same set of time-ordered 
diagrams, but disregarding the thermal character of its external lines: this 
is precisely the set of time-ordered graphs in zero-temperature theory that 
add up to a single, fully covariant, Feynman diagram with the given set of 
external lines. Then the rules above, apart from the product of thermal 
factors (d.2) (which is the same for all the diagrams in the set), provide an 
alternative way for computing this very same zero-temperature Feynman di- 
agram. In fact, it is easy to see that rule (d.l) gives the same contribution, 
regardless of whether we consider some of the thermal legs as internal or all 
of them as external. In our simple example above, for instance, equation (TO) 
corresponds to the regrouping of graphs l.c and l.f. Notice that the thermal 
legs attach in identical manners to the vertices in both graphs. 

Rules for a simpler approach: There is therefore an alternative simpler 
way of computing a given Euclidean finite-temperature graph: 

a. Generate all possible connected graphs that can be obtained from the 
given graph by snipping any number of internal lines, and attach crosses 
to the ends of the snipped lines (the thermal legs). For any given 
thermal pair, one of the legs has to be incoming and the other one 
outgoing, both carrying the same momentum and energy (— iEi, hi). 
Both possible orientations have to be considered, each one generating 
a different diagram. 

b. Compute the graphs in (a) using standard Euclidean Feynman rules as 
in zero-temperature and integrate each graph over the spatial momenta 
of its thermal pairs, with weight factor ni/2Ei for each pair. 
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The given Euclidean graph is equal to the sum of the contributions above 
plus its zero-temperature value. 

Final Comments: The modifications to our rules for the fermionic case 
are rather simple. First, there is the sacred —1 factor for each fermion loop. 
Second, for each fermionic (snipped or not) line of spatial momentum p and 
mass rrif, the numerator gets a factor rrif, where P s = (—siE p , p), s = ±1, 
and where the Euclidean gamma matrices 7^ (given by = —25^) 

must be used. The variable s takes on the value +1 if the fermionic line 
is oriented in the forward time direction and —1 otherwise. Finally, the 
thermal factor associated with a snipped fermionic line is —n(Ek), where 
n(E) = (exp(E /T) + l)" 1 is the Fermi-Dirac thermal occupation number. 

We should point out that our method, although useful in explicit calcu- 
lations, may not be convenient for renormalization group analysis. In T=0 
Feynman diagrams, low order corrections to self-energies and vertices are 
computed only once; the results can then be inserted into higher order skele- 
tons, which makes the renormalization group simple to implement. This 
same principle works for both the Matsubara method (before continuation) 
and the real-time method with doubled components. However, this feature 
may not work here. 

As a final remark, we should point out a possible extension of our method 
to the real-time formalism. If one considers all Euclidean external energies P ; ° 
as imaginary, such that iP® = E\ becomes real (i.e. an analytic continuation), 
includes an ie prescription such that the denominator in the rule d.l becomes 
of the form E inc — E cut + ie, and replaces the factor (— 1)^ by i in rule d.3, 
then at T = the rules precisely reproduce the Green's functions of the real 
time theory. This ie prescription provides a definite analytic continuation of 
the Euclidean result into Minkowski space. We are presently investigating 
whether this continuation is indeed the real-time theory at finite temperature 
as well. 
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Figure captions 

FIG. 1. The lowest order correction to the inverse propagator in the scalar 
g<j) 3 theory, (a) and (b) are the two possible time-ordered graphs at T = 0. 
(c) to (f) are the thermal contributions derived from (a) and (b) by snipping 
the internal line with energy E^: in (c) and (e) the thermal legs are internal] 
in (d) and (f) they are external. We read off E inc as iPo for (c) and (e), 
and iP + E 1 for (d) and (f). Similarly, E cut is E 1 + E 2 for (c), E 2 for (d), 
2iP + E 1 + E 2 for (e), and 2iP + E 2 for (f). 

FIG. 2 A 2-loop contribution in A0 4 . (a) and (b) are the two possible time- 
ordered graphs at T = 0. (c) and (d) are some of the thermal diagrams 
generated from (a) and (b) by snipping one line: in (c) the thermal legs 
are external; in (d) they are internal, (e) and (f) are some of the diagrams 
generated from (a) and (b) by snipping two lines: in (e) both thermal pairs 
are external; in (f) pair 1 is internal and pair 3 external. 

FIG. 3. One-loop contribution to the 3-point function in g<p 3 theory, (a) and 
(b) are two possible time orderings at T = 0. (c) and (d) are some of the 
thermal diagrams generated from (a) and (b) by snipping one line: in (c) the 
thermal legs are external; in (d) they are internal. 

FIG. 4. (a) t-channel and (b) s-channel Feynman diagrams for the T = 
forward scattering process 4>(P°, p) + (p(—iE 1 , k) — > (f>(P°, p) + (f)(—iE 1 , k). 
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